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A natural framework to understand the energy dependence of bulk observables from lower energy 
experiments to the LHC is provided by the Color Glass Condensate, which leads to a "geometrical 
scaling" in terms of an energy dependent saturation scale Qs. The measured charged multiplicity, 
however, seems to grow faster (~ -v/s" i^i nucleus-nucleus collisions than it does for protons 
(~ y/s^''^), violating the expectation from geometric scaling. We argue that this difference between 
pp and AA collisions can be understood from the effect of DGLAP evolution on the value of the 
saturation scale, and is consistent with gluon saturation observations at HERA. 

PACS numbers: 24.85.+p,25.75.-q,12.38.Mh 



I. INTRODUCTION 

The unprecedented high energies of the LHC proton 
and nuclear beams provide us with new experimentals 
tests of QCD dynamics at high energy. On a fundamen- 
tal level we know that also the bulk properties of the 
collision system such as the momentum spectra and cor- 
relations of all produced hadrons must follow from QCD. 
How this happens and to what extent the process can be 
understood in a weak-coupling approximation is still an 
open issue. However, the phenomenological success of the 
Color Glass Condensate (CGC, for reviews see e.g. [l|) 
makes one optimistic that this could in fact be possible. 
In the CGC framework the small x degrees of freedom 
that dominate bulk particle production in hadronic colli- 
sions are described as nonperturbatively strong classical 
color fields. This picture leads naturally to the concept 
of gluon saturation and the saturation scale Qs as the 
dominant transverse momentum scale determining both 
the magnitude and the space and time dependence (i.e. 
the momenta of the gluons) of the small x gluon fields. 

The first LHC observable, in both proton proton and 
nucleus-nucleus collisions, to give us information about 
QCD dynamics at high energy is the charged particle 
multiplicity [2-4]. As we shall discuss in more detail in 
Sec. El the charged multiplicity is, to a very good ap- 
proximation, proportional to Q^. Thus the energy de- 
pendence of the multiplicity is an experimental probe of 
the X dependence of the saturation scale Q^, separately 
for nucleons and nuclei. The simplest model-independent 
way to see this is to realize that gluon saturation turns 
particle production into a one scale problem, with Qs as 
the only dimcnsionful scale apart from the size of the 
system. All bulk quantities, such as the multiplicity or 
the transverse energy, and correlations in the system can 
be understood in this way up to normalization constants 
parametrically of order 1. 

In the CGC framework the energy dependence of the 
saturation scale follows from the JIMWLK [5] rcnornial- 
ization group equation or its mean field, large Nc approx- 
imation, the BK equation |6|. At an intermediate scale. 
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FIG. 1: The charged multiplicity in pp and A A collisions 
estimated from the saturation scale in the IPsat model. Ex- 
perimental datapoints from Refs. [2, @, 0, HO] for pp and 
Refs. 0, 113, m for AA collisions. For details on relating 
Qs to the charged multiplicity see discussion in Sec. FVl 



typically taken as x 0.01, where one would start the 
BK or JIMWLK evolution, the typical nuclear satura- 
tion scale could be estimated as Q^a ~ cA^^^^Q^p with 
c ~ A^^^Rp/R^ a constant somewhat smaller than 1 (for 
a more detailed discussion see e.g. [jI). With a fixed cou- 
pling constant leading order JIMWLK or BK evolution 
would preserve the value of Qsa/Qsp at all energies. In- 
cluding the running of as, however, changes this picture. 
At asymptotically high energies running coupling BK or 
JIMWLK evolution leads to a saturation scale that is in- 
dependent of A Q and therefore grows more slowly for 
nuclei than for nucleons. 

The arguments outlined above would lead to a primary 
gluon multiplicity that would have the same energy de- 
pendence in pp and AA collisions (fixed coupling) or that 
grows more slowly in AA than pp collisions (running cou- 
pling at asymptotical energies). The trend seen in the 
data on the final charged multiplicity is, however, the 
opposite one, posing a puzzle for attempts to understand 
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it on the basis of gluon saturation. We argue in this paper 
that at least a part of the explanation lies in transient 
effects that do not follow directly from BK/JIMWLK 
evolution. We will show that, due to effects of the nu- 
clear geometry and of DGLAP evolution, the growth of 
Qs with energy can actually be faster in nuclei than in 
protons. Thus a single parametrization of the dipole cross 
section can be in good agreement with the basic features 
of the experimental data in both protons and nuclei, as 
shown in Fig. [T] We will first discuss the energy depen- 
dence in different fits to HERA data in Sec. [IT] before 
looking more closely at what happens in two particular 
dipole cross section parametrizations, "IPsat" in Sec. IIIII 
and "bCGC" in Seclvl We will then, in Sec. El discuss 
in more detail the the relation between the initial glu- 
onic and final charged hadron multiplicity in nuclar and 
proton collisions. 



II. ENERGY DEPENDENCE OF THE 
SATURATION SCALE 
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FIG. 2: The saturation scale (adjoint representation) at the 
median impact parameter femod in a lead nucleus and a proton 
in the IPsat and bCGC models, as discussed in the text. Also 
shown are the values A obtained by fitting an exponential 
Q? = ax~^ to the curves. 



Let us now discuss what is known about the value 
of the saturation scale based on fits to HERA data 
There exists by now a large amount of different saturated 
parametrizations of the dipole cross section, mostly fit to 
HERA or in some cases hj to RHIC dAu data. Usu- 
ally the impact parameter dependence is factorized into a 
Gaussian profile multiplying the dimensionless scattering 
amplitude: this is the case in the well known GBW [ij] 
and HM [IJ] parametrizations and the more recent so- 
lutions of the running coupling BK equation (l5| . In 
these cases the fit seems to favor an x-dependence of the 
saturation scale that is faster than the observed energy 
dependence of the charged multiplicity in pp collisions. 
This is the case also for the BK evolution studied in [l^ , 
where the evolution speed is not a fit parameter, but fol- 
lows directly from the evolution equation itself. Thus 
the value A ~ 0.3 obtained from a simple power law 
fit dNch/ dry ^ , to the energy dependence of the 
multiplicity in AA collisions is in good agreement with 
the evolution speed obtained e.g. in the GBW 13] fit 
to HERA data. Although some saturation calculations 
slightly underpredicted the AA multiplicity at the LHC, 
they still correctly predict a stronger growth than in pp 
collisions. The more recent application of running cou- 
pling BK evolution, with a more detailed inclusion of 
the nuclear geometry, reproduces the ALICE multiplic- 
ity data perfectly [161]. 

On the other hand, it appears that in dipole cross sec- 
tion parametrizations where the impact parameter de- 
pendence is not factorized out but included in the sat- 
uration scale itself, fits to HERA data prefer a slower 
increase of Qs with x. This is the case in both of the 
parametrizations, IPsat [13, \M and bCGC fll [11] that 
we shall analyze in more detail in Sees IIIII and IIVI This 
is due to the functional form that intertwines the r and 
b dependence of the dipole cross section in such a way 



that the proton is allowed to grow with energy; leading 
to a growing total DIS cross sections even with a slower 
increase of the saturation scale. This growth is consis- 
tent with the increase of the total pp cross section with 
energy and the t-dependence of diffractive vector meson 
production at HERA. As discussed e.g. in [2^, this 
is the kind of parametrization that one would generally 
prefer on theoretical grounds, since it causes the dipole 
cross section to saturate towards the correct black disk 
limit also for 6^0. However, it is not obvious if the 
particular functional form chosen in these parametriza- 
tions is the correct one. The slower growth of Qs in the 
IPsat and bCGC parametrizations has, in /cT-factorized 
calculations, yielded a good agreement with the charged 
multiplicity in LHC proton-proton collisions l2l|. Also 
the multiplicity distributions and pT-spectra in pp colli- 
sions have recently been analyzed in the /cT-factorization 
approach [22] . 

In conclusion, among the different CGC fits to HERA 
data there are ones that explain the multiplicity in pp 
collisions, and others that give a good description of the 
multiplicity in A A collisions. Neither the pp or A A mul- 
tiplicity data separately is thus an indication against 
gluon saturation, but the apparent failure to describe 
both with the same parametrization is problematic. As 
we shall now see, this is in fact not the case. The IPsat 
parametrization, and also the bCGC (with an additional 
assumption that we will discuss) , can in fact describe the 
^/s dependence of both the pp and AA multiplicities. 



III. NUCLEAR EFFECTS IN EIKONALIZED 
DGLAP 

The IPsat model [13, [3 a modification of the idea 
(see e.g. ^]) of including multiple scatterings and enforc- 
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FIG. 3: The function F{x,r) of Eq. ([2]), proportional to the 
DGLAP-evolved gluon distribution, appearing in the IPsat 
model. 



ing the black disk limit in the DIS cross section by expo- 
nentiating (eikonalizing) a DGLAP-evolved gluon distri- 
bution. In the IPsat parametrization the impact param- 
eter dependence is included in the saturation scale (or 
DGLAP gluon distribution), leading to the dipole cross 
section 



d2bT 



(bT, r, x) = 2 [1 - exp {-r^F{x, r)Tp{h))] . (1) 



Here Tp(b) = exp (— &2/2i3p) /{2'kB) is the impact pa- 
rameter profile function in the proton with Bp = 
4.0 GeV^ and F is proportional to the gluon distribu- 
tion 

On a conceptual level, this formulation resums the mul- 
tiple scatterings off the small- a; gluons by assuming that 
they are independent and thus exponentiate. Strictly 
speaking it is not a "CGC" parametrization in the 
sense that it does include the correlations between the 
small-x gluons in the target that are included in the 
JIMWLK/BK equations. On the other hand it has the 
correct DGLAP behavior at large (small r) and should 
be seen as a a good way to approach the saturation 
regime from the high or large x direction. 

A generalization of Eq. ^ to nuclei (including fluctu- 
ations in the positions of the nucleons) is very straight- 
forward: one replaces the thickness by a sum over A nu- 
cleons as 



= 2 



1-expi -r2F(a:,r)^rp(bT-bT,) 



i=l 



d^br 

This corresponds to treating the interactions with the 
separate nucleons in the nucleus as independent; con- 
sistently with the eikonalization idea. In terms of the 



S'-matrix of the dipole scattering off the target Eq. ([3]) is 
equivalent to 



SA{r, by, a;) = Sp{r, hr - hTi,x). 



(4) 



The positions of the nucleons can then be averaged over 
to yield a Glauber-like averaged dipole cross section, 
written in the large A approximation as 



d^dip 
d^br 



1 — exp 



ATAihr) 



'dip 



(r,x) 



(5) 



where o'dip('"'^) ^'^^ nucleon dipole cross section of 
Eq. (IlJ integrated over the impact parameter by. 

The saturation scale Qs characterizes the qualitative 
change between the dilute color transparency region r — > 
and the black disk limit -grg^ — > 2 at large r. We 
shall take here a mo del- independent definition of Qs as 
the solution of 



dcTdip 

d^br 



(x, r - l/Qs{x, bT)) = 2(1 - e-i/4). (6) 



The saturation scales defined as Eq. © for a proton and 
a lead nucleus in the IPsat parametrization are plotted in 
Fig. [21 multiplied by the color factor Ca/Cf appropriate 
for gluon production in pp or AA collisions. We are tak- 
ing an effective value representing the average over the 
transverse plane by taking the saturation scale at 6mcd, 
the median impact parameter of the total DIS cross sec- 
tion (the value such that half the cross section comes 
from h < ferncd)- The difference between protons and nu- 
clei is striking: the energy dependence is • 
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nuclei and 



-0.20 



for 

for protons, in a manner which 
immediately evokes the y^-dependence of the charged 
multiplicity. 

The reason for this difference lies in the behavior of the 
DGLAP-evolved gluon distribution, whose x-dependence 
gets steeper at higher Q and, because Qs grows with A, at 
higher A. This feature was mentioned already in Ref. [7|, 
although the discussion there is formulated in terms of 
the A-dependence at fixed x in stead of the x-dependence 
at fixed A. The function F{x, r) of Eq. ^ is shown in 
Fig. [3] The gluon distribution at the initial scale 
in the IPsat model has a very mild x-dependence. The 
DGLAP evolution then drives the distribution to become 
much steeper at higher scales which, because of the A^^^ 
enhancement of Qg in nuclei, define the saturation region 
in a nucleus. 

Let us try to estimate the saturation scales explicitly 
to illustrate how this happens. From Eqs. (HJ, ^ and 
([6]) we can estimate the proton saturation scale at the 
center of the proton as 



{b- 



lF(x,r 



l/Qsp{b = 0)) 



2ttB„ 



(7) 



Replacing, for purposes of illustration, the Woods-Saxon 
profile by a theta function TA{b) ~ d{RA — b)/{TTR\) we 
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get for the saturation scale in a nucleus 



l/0sA,x)/2 



ttRI 



1 

4' 



(8) 



Now we know that Qsa > Qsp and therefore r = 1/Qsa 
is in the dilute region for the proton dipole cross section. 
We can therefore take in Eq. ([5]) the small r approxima- 
tion 

^dipl'^ = x)/2 « -^Fix, r = 1/Q,a) (9) 

and thus 

Q^.(*^0)44a££zWMfci«i^ (10) 

Equations ([7]) and pHl) are still implicit equations that 
must be solved to obtain the saturation scales. It is, 
however, easy to see that because QsA > Qsp, the gluon 
distribution on the r.h.s. of Eqs. Q and (fTO|). and conse- 
quently the saturation scale on the l.h.s., evolves more 
rapidly with energy. We emphasize that this discus- 
sion is just an illustration of the origin of the different 
x- dependences in protons and nuclei, and the values in 
Fig. [2] are obtained from the full expressions. 

Finally, relating these saturation scales to the charged 
multiplicities as discussed in Sec. |V] results in Fig. [T] 
There are two curves for protons and nuclei, differing 
by whether one keeps the coupling as fixed at 0.33 or 
whether one allows it to run as as(Qs). We see that the 
agreement with the experimental data is extremely good, 
considering the simplicity of this approach. 



IV. SATURATION SCALE FOR INDEPENDENT 
BK-EVOLVED NUCLEONS 

Whether the impact parameter profile is smooth (as in 
e.g. [l^) or fluctuating [Ig], running coupling BK evolu- 
tion leads to a slower increase of Qs for a higher initial 
value, i.e. for nuclei. Let us here consider an approxi- 
mation that the individual nucleons evolve according to 
BK, and are then combined into a nucleus using the as- 
sumption of independent scattering Eq. ^ . Parametri- 
cally this ansatz could perhaps be justified at most in a 
moderate x regime where evolution does not yet happen 
coherently over the whole nucleus, and would certainly 
not be valid for asymptotically high energies. It is, how- 
ever, interesting to see how a faster energy dependence 
in nuclei can arise also in this scenario. 

To be more precise, the dipole cross section for a proton 
in the bCGC parametrization [l^, [l^ is: 



d^b 



for rQ'^ < 2 



The saturation scales Qs and Q'^ are conceptually the 
same quantity and their numerical values are of the 
same order, but we differentiate between them in order 
to maintain our model independent definition of Qs in 
Eq. (|6]). The coefficients A and B in the can be de- 
termined uniquely from the condition that -grg^ and 
its first derivative with respect to rQ'^ are continuous 
across rQ'^ = 2. Here Y = ln(l/x) is the rapidity, while 
7s = 0.628 and k = 9.9 (which quantifies the geometric 
scaling violations in Eq. pTjl ) are obtained from leading 
logarithmic BFKL dynamics [2j| . The impact parameter 
dependence of the proton saturation scale is introduced 
into the bCGC model in the form 

Q'.ix, h) = (^) ^ [exp (-6V2i?cGc)] ^ GeV . (12) 

The parameters A, xq, Aq and Bqgc s^re fit to the data, 
with the fit resulting in A = 0.159, xq = 5.95 • 10"'', 
Afo = 0.417 and BcGC = 5.5 GeV"^ 

We now use this parametrization for nuclei by assum- 
ing that the scatterings off the nucleons are independent, 
i.e. assuming Eq. (jlj which leads to the Glauber form for 
the avergare gluon distribution ([3]). The resulting satu- 
ration scales, at the median impact parameter b^cd are 
plotted in Fig. [2] for protons and nuclei. Due to the non- 
trivial functional form of the parametrization, the evolu- 
tion speed for protons turn out to be 0.18, slightly larger 
than the parameter A in the parametrization, but still 
slower than in 6-independent fits or in BK evolution. The 
nuclear saturation scale grows faster, as Q^a ~ x~^'^^. 

Let us try to understand this difference analytically in 
the similar way as in Sec. IIIII It is easier here to use the 
definition of the saturation scale Qi, that appears in the 
parametrization itself. Thus we have for the proton 



-{r 



d^br 
and in the nucleus 



d^hr 



2 



2/Q'x)^2Uo 



(13) 



1 — exp 



ATA{b) 



'dip 



(14) 



We again replace the Woods-Saxon distribution by a 
theta function Tyi(6) « 9{Ra — b)/{'KR\). Because at 
the nucleus saturation scale one is in the dilute regime 
for the proton we can now approximate Eq. ()14[) by 



Mo 



ATA{b) 



{r^llQ',^,x) 



(15) 



We now assume that the integral over the impact parame- 
ter in the proton approximately factorizes into a constant 
CTo ~ 27ri?cGC to get 



2-2exp(-^ln2(Srg^)) for rQ^ > 2 . '^dipl'^ = VQsA, 2;) « 2(7oA/'o 

(11) 



Q'sA 



2 7. + 7 



(16) 



leaving us with 
ttRI 



exp 



7s + 



1 



2kXY 



A A 



where we have denoted 

0' 



sA 



i sp 



„A/2 



(17) 



(18) 



Since the l.h.s of Eq. (fTT]) is independent of energy, it is 
obvious that A must grow with the energy (i.e. with Y). 
Differentiating with respect to Y gives 



A'(y) = 



A2 



2kXY^j, 



A ^ 



(19) 



Assuming A w \nA^/^ w 1.8 and Y = In(lOOO) w 7 
we get A' « 0.03. In terms of the saturation scales 
this means that ~ Q' x~^-°^, which explains most 
of the effect seen in Fig. ^ The interpretation of this 
result is in fact the same as in the IPsat case. In the 
bCGC parametrization there is a logarithmic term in 
the exponent that violates geometric scaling. At smaller 
Y, i.e. larger x, the effective anomalous dimension 
7eff = 7s + ln(2/(r(3g))/(KAF) is larger, i.e. closer to 
1; thus the Q^-dependence of the integrated gluon dis- 
tribution is close to a logarithm. At large Y or small x 
one recovers the anomalous dimension 7^ , which leads to 
a much faster increase of the integrated gluon distribu- 
tion with Q^. This is precisely the scenario that lead to 
a faster growth of QsA in the case of the IPsat model. 
Another way to see this is to rewrite Eq. (fT7|) as 



sp 



f Aao_ 



/3 



(20) 



At smaller x is larger, thus the nuclear enhance- 

ment of Qs is larger. 

Again, using the procedure described in Sec. IVl results 
in the estimates for the charged multiplicity shown in 
Fig. 131 While the agreement with experiental data is not 
as good as with the IPsat parametrization, the general 
trend of a faster increase in AA than in pp is still seen. 



V. RELATION BETWEEN AND iVch 

In any CGC calculation of gluon production in a col- 
lision of two hadronic objects the initial gluon multiplic- 
ity depends on the saturation scale parametrically in the 
same way. The theory of the CGC is based on weak cou- 
pling calculations, and for the consistency of the frame- 
work one assumes that as(<9s) 1- This means that Qs 
is a semihard scale and we can assume that parametri- 
cally Qs ^ Aqcd- In the CGC, the saturation scale also 
defines the correlation length of the system in the trans- 
verse plane, ~ 1/Qs- In the limit when a weak coupling 
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FIG. 4: Energy dependence of the charged multiplicity esti- 
mated using the bCGC saturation scale, extended to nuclei 
as described in Sec. I IVl 



CGC calculation is justified in the first place, the corre- 
lation length is smaller than the size of the interaction 
region, l/Q^ ^ cr. Thus particle production happens 
locally in independent domains of size ~ 1/Qs in the 
transverse plane. This picture leads to a gluon multiplic- 
ity that can be written as a local observable, where for 
dimensional reasons the number of gluons per unit area 
is proportional to the local Q^: 



dNi 



init.p 



(21) 



Here, following [25|,|26|, we have introduced the "gluon 
liberation coefficient" c; a nonperturbative dimensionless 
constant that is parametrically of order 1, but depends on 
the detailed spectrum of the produced gluons. Its value 
in the MV model [111 has been determined using Classical 
Yang-Mills simulations '^s', '29'] to be c « 1 . 1 (see [30| for 
a discussion of the CYM results parametrized in terms of 
c given here). As discussed e.g. in Ref. [3l| the value of c 
remains very close to this value across a range of different 
models for the color charge distribution (or, equivalently, 
for the dipole cross section). There is also an analytical 
calculation |26| of the liberation coefficient with the result 
CR:! 21n2 « 1.4. 

Phenomenological studies are often done using various 
/cT-factorized approximations to compute the gluon spec- 
trum. Although one can derive a fcT-factorized formula 
for the dilute-dilute "pp" and dilute-dense "pA" cases, 
/cT-factorization yields the wrong gluon s pec trum for the 
case of dense-dense or "AA" scattering jS^]- Neverthe- 
less, since the only dimensionful scale in the problem is 
still Qs, also the result from /cT-factorization can still be 
parametrized as (HH), although the value of the coeffient 
c is incorrect. 

In proton proton collisions, the proportionality be- 
tween the initial gluon and final charged hadron mul- 
tiplicities is based on the phenomenological success of 



6 



local parton-hadron duality (LPHD) [33*1 . The working 
assumption is that the final charged multiplicity is pro- 
portional to the initial partonic one. The constant of 
proportionality is a property of the independent frag- 
mentation of partons into hadrons and should thus be 
basically independent of collision energy or centrality. In 
order to explain the slower growth of the charged mul- 
tiplicity in pp-coUisions by a modification of the LPHD 
hypothesis one would have to argue that as the energy 
increases, the number of charged hadrons produced from 
a primary parton decreases, which does not seem likely. 
Also the transverse area of the interaction region should, 
if anything, increase with energy. Therefore it seems that 
the proton saturation scale should grow with energy at 
most as the observed charged multiplicity, or even more 
slowly. 

On the side of nucleus-nucleus collisions there is a 
wider range of plausible modifications to the relation be- 
tween initial gluonic and observed hadronic multiplici- 
ties. One possible starting point is the LPHD assump- 
tion as in pp collisions. The fragmentation could then 
be argued to lead to a larger ratio of the final charged 
hadron to the initial gluon multiplicity, e.g. due to the 
larger initial gluon momentum {px) Qs] this is the 
argument in e.g. jsj]. One potential problem for this 
approach comes from considering the first moment of 
the particle spectrum, i.e. the total transverse energy. 
For the initial gluons the typical initial transverse mo- 
mentum is {pt) ^ Qs ^ \J ( dA^/ d7;)/5i. A fragmen- 
tation process that leads to a larger final charged mul- 
tiplicity for AA than pp collisions should thus lead to 
(pt)/\/ ( dA^/ Ar])IS^ that decreases towards central col- 
lisions and towards higher energies. 

To be more explicit, let us assume that in the LPHD 
picture one gluon produces n charged particles after frag- 
mentation. To account for the faster growth of the multi- 
plicity with yfs in A A collisions, one would need ua > np. 
We now have for the initial gluons 



{PT)g ~ Qs 

Si_ d^ ^ 
and for the final charged particles 

(Pt) ch ^ Qs/n 
1 dA^ch „2 



(22) 
(23) 

(24) 
(25) 



since transverse momentum must be conserved during 
fragmentation. Now, if iia > np, the scaled mean px 



(pt) 



ch 



1 



1 djVdi ny'n 



(26) 



should be smaller in central AA collisions than for pro- 
tons. This is indeed seen in the RHIC data [lO| . However, 
if the increase in n is due the larger (p-r) ~ Qs the ratio 



(1261) should also decrease with increasing collision energy. 
Between = 62.4 GeV and ^/s = 200 GeV no such de- 
crease is seen at RHIC [l^] , while no firm conclusions for 
higher energies can yet be made from the LHC data. 

The LPHD scenario is more or less based on neglect- 
ing all collective effects even in A A collisions, and thus 
assuming that no quark gluon plasma is formed. Am- 
ple experimental evidence points to the contrary. In the 
extreme case of strong interactions among the produced 
gluons, they thermalize into an isotropic plasma, which 
then expands in the transverse, and, more importantly, 
the longitudinal direction according to (nearly) ideal hy- 
drodynamical equations of motion. In the ideal hydro- 
dynamical case the entropy, and thus multiplicity, of the 
particles stays constant during the evolution. This leads 
to the picture where the final (total) multiplicity is, not 
only proportional, but nearly equal to the initial glu- 
onic one. During the hydrodynamical expansion of a lo- 
cally isotropic system the mean transverse momentum 
(or energy per particle) decreases by a large amount due 
to p dV work done pushing the expanding plasma down 
the beampipe. Radial fiow developing during the evolu- 
tion also boosts the transverse momenta of the particles 
compared to pp collisions, consistently with the observed 
increase with centrality of (pt)- The drawback of this 
scenario is that, in spite of much work, we do not have 
a quantitative theoretical understanding of the thermal- 
ization process. This fast thermalization argument is, 
however, strangthened by the results of explicit calcula- 
tions of the initial gluon multiplicity in the CGC, which 
yield an initial gluon multiplicity that is close to the final 
total multiplicity. The initial transverse energy, on the 
other hand, is larger than the observed one by a large fac- 
tor, which would be consistent with the hydrodynamical 
picture. 

The true physical situation is most likely to lie some- 
where between the two extreme scenarios of LPHD and 
ideal hydrodynamics, with some entropy production and 
thus increase in the multiplicity during the spacetime 
evolution of the plasma. However, based on this discus- 
sion it seems unlikely that final state effects would solve 
the problem of a faster growth of the multiplicity with 
^/s in A A collisions than in pp. 

Let us now use the assumption of fast thermalization 
and ideal hydrodynamical expansion to obtain quanti- 
tative estimates for the charged multiplicity. For central 
heavy ion collisions we take the following simple estimate 
for the final charged multiplicity multiplicity: 



2 dA^ch 2 2 dN„ 



A, 



part 



drj 



3 A, 



part 



2 CfQK^) 2S^ 

3 2TT^as Apart 



Here the factor 2/3 accounts for the fraction of charged 
particles of the total multiplicity (to a first approximation 
Tr^jVr"). We take here c — 1.1 as discussed above. The 
typical transverse area S± per participant pair is taken 
as the value estimated in central gold-gold collisions by 
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STAR m as 



2Sj_ 154 fm^ 



part 



0.5 X 350 



(28) 



The saturation scale as extracted from fits to DIS data 
depends on the momentum fraction x. In hadronic or 
heavy ion colhsions the corresponding variable is a ra- 
tio of the transverse momentum to the collision energy. 
The X and Qs{x) corresponding to each collision energy 
Qs{^/s) are solved from the relation 



Qsjx) 



(29) 



We emphasize that, unlike in typical fc^-factorized cal- 
culations, there is no arbitrary normalization factor to 
adjust here. Once the saturation scale is known, it de- 
termines both the normalization and the shape of the pt 
spectrum of the produced gluons. The mean transverse 
momentum of the initial gluons is (pt) depends slightly 
more on the precise kx dependence of the dipole cross 
section than the multiplicity. For the case of the MV 
model the spectrum is relatively hard, with (px) ~ 1.3Qs 
(following Refs. [H, H^l). For RHIC this corresponds 
to (pt) ~ 1.5 GeV for the intial gluons. With rapid 
thermalization and (nearly) ideal (nearly) boost invari- 
ant hydrodynamical evolution this is then reduced (see 
e.g. [23]) by a factor 3 to match the final observed trans- 
verse energy of around 0.5 GeV per particle. At the 
2.75 A TeV collision energy at the LHC a similar esti- 
mate yields (pr) ~ 2.0 GeV for the gluons in the initial 
state which would be similarly reduced by hydrodynam- 
ical evolution. 

Now this can be contrasted with the estimate based on 
kx factorization [l6j . where, after adjusting the normal- 
ization constant to reproduce the RHIC multiplicity, one 
obtains at the LHC a transverse energy^ 14 TeV and a 
total multiplicity « (3/2)1600 — 2400, i.e. a transverse 
energy of 6 GeV per particle. Although a value for the 
transverse energy has not yet been released by the LHC 
experiments, based on the spectra published in Ref. [1^ 
it seems safe to say that this would require a reduction of 
at least a factor of 6 between the initial gluons and the 
final state particles. This failure to calculate both the 
initial multiplicity and transverse energy without adjust- 
ing both by an arbitrary normalization constant follows 
from the incorrect gluon spectrum in fcT-factorization. 

For protons we simply replace Eq. (j27|) with an ex- 
pression that leaves out the scaling by the number of 
participant pairs 



djVch 
dr/ 



2 CfQUx) ^ 
3^^ 2n^as 



(30) 



Here, adjusting the normalization to data as always with 
LPHD, we take as the transverse area by a constant S± — 
20 mb for fixed as and S±_ = 24 mb for running as- 
This includes the conversion from gluons to final hadrons. 
This summarizes the simple procedure used to arrive at 
the multiplicity estimates in Figs. [T] and |31 



VI. CONCLUSIONS 

The experimental data seems to be hinting that for 
bulk particle production at LHC energies one is not yet 
far enough the asymptotic high energy, Qs ^ Aqcd, 
regime for gluon saturation to work perfectly without 
additional finite -y/s corrections. What these corrections 
are remains still somewhat an open issue. One possibil- 
ity is that soft confinement scale physics remains to be 
dominant in pp collisions; in the weak coupling frame- 
work pursued in this paper we have not been able to ana- 
lyze this option. The other possibility, more encouraging 
in terms of prospects for first principles understanding, 
is that the nuclear saturation scale is large enough at 
the LHC to be sensitive to the large effects such as 
DGLAP evolution. Note that also in the BK equation 
the glowth of the saturation scale with energy is slower 
in the preasymptotic regime close to the initial condition. 
In phenomenological applications [l^, 113] this preasym- 
totic slower growth is essential for agreement with ex- 
perimental data. Thus the energy dependence at LHC 
energies is to a large degree a consequence of the trans- 
verse momentum dependence in the initial condition, not 
only of the evolution itself. 

We have in this note shown that the different energy 
dependences of the charged particle multiplicities in pp 
and AA collisions can be understood in the framework of 
gluon saturation. There is a well-tested and motivated 
impact parameter-dependent parametrization, based on 
an eikonalized, DGLAP-evolved gluon distribution, that 
very accurately describes both pp and AA multiplicities. 
The difference between protons and nuclei comes, in this 
parametrization, from effects of DGLAP evolution on the 
slope of the gluon distribution. These effects are not 
present in the pure JIMWLK/BK evolution formalism, 
whether at fixed or running coupling. Fully understand- 
ing the physics at play here requires incorporating higher 
transverse momentum physics not only into the initial 
condition, but also into nonlinear evolution itself. Some 
steps this direction have already been taken (38| , but fur- 
ther investigation is needed to fully understand the con- 
sequences for bulk particle production at the LHC. 
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